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Abstract 

We construct a statistical model of crystal melting to count BPS bound states of DO and D2 branes 
on a single D6 brane wrapping an arbitrary toric Calabi-Yau threefold. The three-dimensional 
crystalline structure is determined by the quiver diagram and the brane tiling which characterize 
the low energy effective theory of D branes. The crystal is composed of atoms of different colors, 
each of which corresponds to a node of the quiver diagram, and the chemical bond is dictated by the 
arrows of the quiver diagram. BPS states are constructed by removing atoms from the crystal. This 
generalizes the earlier results on the BPS state counting to an arbitrary non-compact toric Calabi- 
Yau manifold. We point out that a proper understanding of the relation between the topological 
string theory and the crystal melting involves the wall crossing in the Donaldson-Thomas theory. 



1 Introduction 



In type IIA superstring theory, supersymmetric bound states of D branes wrapping holomorphic 
cycles on a Calabi-Yau manifold give rise to BPS particles in four dimensions. In the past few 
years, remarkable connections have been found between the counting of such bound states and the 
topological string theory: 

(1) When the D brane charges are such that bound states become large black holes with smooth 
event horizons, the OSV conjecture [1] states that the generating function Zbh of a suitable index 
for black hole microstates is equal to the absolute value squared of the topological string partition 
function Z^^pi 

Zbk = \Ztop\^ , (1-1) 
to all orders in the string coupling expansion. 

(2) When there is a single D6 brane with DO and D2 branes bound on it, it has been proposed [2] that 
the bound states are counted by the Donaldson-Thomas invariants [3,4] of the moduli space of ideal 
sheaves on the D6 brane. For a non-compact toric Calabi-Yau manifold, the Donaldson-Thomas 
invariants are related to the topological string partition function [2, 5, 6] using the topological 
vertex construction [7]. Recently the connection between the topological string theory and the 
Donaldson-Thomas theory for toric Calabi-Yau manifolds was proven mathematically in [8] . Given 
the conjectural relation between the counting of D brane bound states and the Donaldson-Thomas 
theory, it is natural to expect the relation, 

Zbu — Ztop- (1-2) 

The purpose of this paper is to understand the case (2) better. We start with the left-hand 
side of the relation, namely the counting of BPS states. Recently, the non-commutative version 
of the Donaldson-Thomas theory is formulated by Szendroi [9] for the conifold and by Mozgovoy 
and Reineke [10] for general toric Calabi-Yau manifolds^ In this paper, we will establish a di- 
rect connection between the non-commutative Donaldson-Thomas theory and the counting of BPS 
bound states of DO and D2 branes on a single D6 brane. Using this correspondence, we will find a 
statistical model of crystal melting which counts the BPS states. 

The crystal melting description has been found earlier in the topological string theory on the 
right-hand side of (|1.2p . It was shown in [2,5] that the topological string partition function on 
C'^, the simplest toric Calabi-Yau manifold, and the topological vertex can be expressed as sums of 
three-dimensional Young diagrams, which can be regarded as complements of molten crystals with 

^See [11—16] for further developments. 
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the cubic lattice structure^ Since the topological vertex can be used to compute the topological 
string partition function for a general non-compact toric Calabi-Yau manifold, it is natural to expect 
that a crystal melting description exists for any such manifold. To our knowledge, however, this 
idea has not been made explicit. The crystal melting model defined in this paper appears to be 
different from the one suggested by the topological vertex construction. 

The low energy effective theory of DO and D2 branes bound on a single D6 brane is a one- 
dimensional supersymmetric gauge theory, which is a dimensional reduction of an A/" = 1 gauge 
theory in four dimensions. The field content of the gauge theory is encoded in a quiver diagram 
and the superpotential can be found by the brane tiling [25-28] cl From these gauge theory data, 
we define a crystalline structure in three dimensions. The crystal is composed of atoms of different 
colors, each of which corresponds to a node of the quiver diagram and carries a particular combi- 
nation of DO and D2 charges. The chemical bond is dictated by the arrows of the quiver diagram. 
There is a special crystal configuration, whose exterior shape lines up with the toric diagram of the 
Calabi-Yau manifold. Such a crystal corresponds to a single D6 brane with no DO and D2 charges. 
We define a rule to remove atoms from the crystal, which basically says that the crystal melts from 
its peak. By using the non-commutative Donaldson-Thomas theory [10], we show that there is a 
one-to-one correspondence between molten crystal configurations and BPS bound states carrying 
non-zero DO and D2 charges. The statistical model of crystal melting computes the index of D 
brane bound states. 

The number of BPS states depends on the choice of the stability condition, and the BPS count- 
ings for different stability conditions are related to each other by the wall crossing formulae. In this 
paper, we find that, under a certain stability condition, BPS bound states of D branes are counted 
by the non-commutative Donaldson-Thomas theory. We can use the wall crossing formulae recently 
derived in [13, 14] to relate this result to the commutative Donaldson- Thomas theory. Since the 
topological string theory is equivalent to the commutative Donaldson- Thomas theory for a gen- 
eral toric Calabi-Yau manifold [8], the relation (|1.2p is indeed true for some choice of the stability 
condition, as expected in [2]. 

In general, the topological string partition function and the partition function of the crystal 
melting model are not identical, but their relation involves the wall crossing, 

^crystal melting ~ -^top (modulo waU CrOSSiugs). (1.3) 

This does not contradict with the result in [2,5] since there is no wall crossing phenomenon for C^. 
In general, however, a proper understanding of the relation between the topological string theory 
and the crystal melting requires that we take the wall crossing phenomena into account. 

^ See [17-24] for further developments. 

^See [29,30] for reviews of the quiver gauge theory and the brane tihng method. 
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In section 2, wc will summarize the computation of D branc bound states from the gauge theory 
perspective. In section 3, we will discuss how this is related to the recent mathematical results on 
the non-commutative Donaldson-Thomas invariants. In section 4, we will formulate the statistical 
model of crystal melting for a general toric Calabi-Yau manifold. The final section is devoted to 
summary of our result and discussion on the wall crossing phenomena. The Appendix explains the 
equivalence of a configuration of molten crystal with a perfect matching of the bipartite graph. 

2 Quiver Quantum Mechanics 

In the classic paper by Douglas and Moore [31], it was shown that the low energy effective theories 
of D branes on some orbifolds are described by gauge theories associated to quiver diagrams. Sub- 
sequently, this result has been generalized to an arbitrary non-compact toric Calabi-Yau threefold. 
A toric Calabi-Yau threefold Xa is a fiber bundle of x M over M?, where the fibers are special 
Lagrangian submanifolds. The toric diagram A tells us where and how the fiber degenerates. For 
a given X^. and a set of DO and D2 branes on X^, the following procedure determines the field 
content and superpotential of the gauge theory on the branes. We will add a single D6 brane to 
the system later in this section. 

2.1 Quiver Diagram and Field Content 

The low energy gauge theory is a one-dimensional theory given by dimensional reduction of an 
jV = 1 supersymmetric gauge theory in four dimensions. The field content of the theory is encoded 
in a quiver diagram, which is determined from the toric data and the set of D branes, as described 
in the following. A quiver diagram Q = (Qg, Qi) consists of a set Qq of nodes, with a rank TVj > 
associated to each node i G Qo, and a set Q\ of arrows connecting the nodes. The corresponding 
gauge theory has a vector multiplet of gauge group U{Ni) at each node i. There is also a chiral 
multiplet in the bifundamental representation associated to each arrow connecting a pair of nodes. 

In the following, we will explain how to identify the quiver diagram. The reader may want 
to consult Figure 1, which describes the procedure for the Suspended Pinched Point singularity, 
which is a Calabi-Yau manifold defined by the toric diagram in Figure l-(a) or equivalently by the 
equation, 

xy = zw^, (2.1) 

in C^. 

To identify the quiver diagram, we take T-dual of the toric Calabi-Yau manifold along the 
fibers [26,32]. The fibers degenerate at loci specified by the toric diagram A, and the T-duality 
replaces the singular fibers by NS5 branes [33]. Some of these NS5 branes divide into domains as 
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Figure 1: (a) The toric diagram for the Suspended Pinched Point singularity, (b) The configuration 
of D2 and NS5 branes after the T-duahty on T^. The green exterior Hues are periodically identified. 
The red lines representing NS5 branes separate the fundamental domain into several domains. The 
T-dual of DO branes wrap the entire fundamental domain, and fractional D2 branes are suspended 
between the red lines. The white domains contain D2 branes only. In each shaded domain, there is 
an additional NS5 brane. There are two types of shades depending of the NS5 brane orientation. 
The white domains are connected by arrows through the vertices, and the directions of the arrows 
are determined by the orientation of the NS5 branes. (c) The quiver diagram obtained by replacing 
the white domains of (b) by the nodes. 



shown in the red lines in Figure l-(b) [30,34-36]. The DO branes become D2 branes wrapping the 
whole T^. The original D2 are still D2 branes after the T-duality, but each of them is in a particular 
domain of T'^ suspended between NS5 branes. In addition, there are some domains that contain 
NS5 branes stretched two-dimensionally in parallel with D2 branes 3 Let us denote the domains 
without NS5 branes by i G Qo and the domains with NS5 branes by a G /. In Figure l-(b), the 
Qo-typc domains are shown in white, and the /-type domains are shown with shade. There are two 
types of shades, corresponding to two different orientations of NS5 branes. This distinction will 
become relevant when we discuss the superpotential. 

The Qo-type domains are identified with nodes of the quiver diagram since open strings ending on 
them can contain massless excitations. The rank Ni of the node i 6 Qo is the number of D2 branes 
in the corresponding domain. On the other hand, /-type domains give rise to the superpotential 
constraints as we shall see below. Though two domains i,j G Qo never share an edge, they can 
touch each other at a vertex. In that case, open strings going between i and j contain massless 
modes. We draw an arrow from i ^ j or i <— j depending on the orientation of the massless open 
string modes, which is determined by the orientation of NS5 branes. Note that the quiver gauge 
theory we consider in this paper are in general chiral. This completes the specification of the quiver 
diagram. 

As another example, the quiver diagram for the conifold geometry has two nodes connected by 
two sets of arrows in both directions. The ranks of the gauge groups are hq and uq — n2, where no 



* The NS5 branes are also filling the four-dimensional spacetime M^''^ while the D2 branes are localized along a 
timelike path in four dimensions. 
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and n2 are the numbers of DO and D2 branes. The gauge theory is a dimensional reduction of the 
Klebanov-Wittcn theory [37] when 712 = and the Klebanov-Strassler theory [38] when 712 > 0. 

2.2 Superpotential and Brane Tiling 

Each domain a G / containing an NS5 brane is surrounded by domains ii,?2, S Qo without 
NS5 branes, as in Figure l-(b). By studying the geometry T-dual to Xa in more detail, one finds 
that the domain is contractible. Since open strings can end on the domains ii, 12, im the domain 
a can give rise to worldsheet instanton corrections to the superpotential. 

This fact, combined with the requirement that the moduli space of the gauge theory agrees with 
the geometric expectation for D branes on Xa, determines the superpotential. Depending on the 
NS5 brane orientation, the /-type domains are further classified into two types, /+ and and 
thus the regions of torus is divided into three types Qo, 1+ smd Such a brane configuration, 
or a classification of regions of T^, is called the brane tiling^ In Figure l-(b), the brane tiling is 
shown by the two different shades. The superpotential W is then given by 

^ = E Tr (li A- .<;>,+ ) - E Tr (n A- l ' ^'''^ 

ael+ \q=l I ael- \q=l ) 

where the domain a E I± are surrounded by the arrows i^^j. — ^ 4°± ^ ■ • • ^ ± ~* ^i"!" 

For each arrow —> ± (1 < 9 < ), the corresponding bifundamental field is denoted 

by A. (a) .(a) . This formula is tested in many examples. In particular, it has been shown that 

\,±'%+i,± 

the formula reproduces the toric Calabi-Yau manifold as the moduli space of the quiver gauge 
theory [40]. 

In the literature of brane tiling, bipartite graphs are often used in place of brane configurations 
as in Figure l-(b). A bipartite graph is a graph consisting of vertices colored either black or white 
and edges connecting black and white vertices. Since bipartite graphs will also play roles in the 
following sections, it would be useful to explain how it is related to our story so far. For a given 
brane configuration, we can draw a bipartite graph on as follows. In each domain in /+ (/-), 
place a white (black) vertex. Draw a line connecting a white vertex in a domain i £ /+ and a black 
vertex in a neighboring domain j G /-. The resulting graph F is bipartite. See Figure 2 for the 
comparison of the brane configuration and the bipartite graph in the case of the Suspended Pinched 
Point singularity. We can turn this into a form that is more commonly found in the literature, for 
example in [26], by choosing a different fundamental region as in Figure 3. 

^In the literature the word brane tiUng refers to the bipartite graph explained below. Here the word brane tiling 
refers to a brane configuration as shown Figure l-(b). Such a graph is called the fivebrane diagram in [39]. 
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brane configuration 



bipartite graph 



Figure 2: The correspondence between the brane configuration on and the bipartite graph. 
The white (black) vertex of the bipartite graph corresponds to the region /+ (/_) in hght (dark) 
shade. The edge of the bipartite graph corresponds to an intersection of /_ and in /+ . From this 
construction, it automatically follows that the graph so obtained is bipartite. 




Figure 3: By choosing a different fundamental region of T^, we find a bipartite graph which is more 
commonly found in the literature. 
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2.3 D-term Constraints and the Moduli Space 

The F-term constraints are given by derivatives of the superpotential, which can be determined as 
in the above. The moduh space of solutions to the D-term constraints is then described by a set of 
gauge invariant observables divided out by the complexified gauge group Gc [41]. The theorem by 
King [42] states that an orbit of Gc contains a solution to the D-term conditions if and only if we 
start with a point that satisfies the ^-stability, a condition defined in the next section. Thus, we 
can think of the moduli space as a set of solutions to the F-term constraints obeying the ^-stability 
condition, modulo the action of Gc- 

2.4 Adding a Single D6 Brane 

To make contact with the Donaldson-Thomas theory, we need to include one D6 brane. Since the 
D6 brane fills the entire Calabi-Yau manifold, which is non-compact, it behaves as a flavor brane. 
In the low energy limit, the open string between the D6 brane and another D brane gives rise to one 
chiral niultiplet in the fundamental representation for the D brane on the other end. The D6 brane 
then enlarges the quiver diagram by one node and one arrow from the new node. To understand 
why we only get one arrow from the D6 brane, let us take T-duality along the fiber again. The 
D6 brane is mapped into a D4 brane which is a point in some region in . This means that we only 
have one new arrow from the new node corresponding to the D6 brane to the node corresponding 
to the D2 branes in the region. See [23,43] for related discussion in the literature. 

3 Non-commutative Donaldson-Thomas Theory 

In the previous section, we discussed how to construct the moduli space of solutions to the F-term 
and D-term constraints in the quiver gauge theory corresponding to a toric Calabi-Yau manifold 
Xa with a set of D0/D2 branes and a single D6 brane. In this section, we will review and interpret 
the mathematical formulation of the non-commutative Donaldson-Thomas invariant in [9, 10] for 
Xa. We find that it is identical to the Euler number of the gauge theory moduli space. 

3.1 Path Algebra and its Module 

For the purpose of this paper, modules are the same as representations. Consider a set of all open 
paths on the quiver diagram Q = {Qo,Qi)- By introducing a product as an operation to join a 
head of a path to a tail of another (the product is supposed to vanish if the head and the tail 
do not match on the same node) and by allowing formal sums of paths, the set of open oriented 
paths can be made into an algebra CQ called the path algebra. We would like to point out that 
there is a one-to-one correspondence between a representation of the path algebra and a classical 
configuration of bifundamental fields of the quiver gauge theory. Suppose there is a representation 
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M of the path algebra. For each node i G Qo, there is a trivial path e-i of zero length that begins 
and ends at i. It is a projection, (ci)^ = e^. Since every path starts at some node i and ends at some 
node j, the sum et acts as the identity on the path algebra. Therefore, M = (BieQgMi, where 
Mi — CiM. Let us write Ni — dim A/;. For each path from i to j, one can assign a map from Mi 
to Mj. In particular, there is an Ni x Nj matrix for each arrow i ^ j G Qi of the quiver diagram. 
By identifying this matrix as the bifundamental field associated to the arrow i — + j, we obtain 
a classical configuration of bifundamental fields with the gauge group U{Ni) at the node i. By 
reversing the process, we can construct a representation of the path algebra for each configuration 
of the bifundamental fields. 

3.2 F-term Constraints and Factor Algebra A 

Let us turn to the F-tcrm constraints. Since the bifundamental fields of the quiver gauge theory is 
a representation of the path algebra, the F-term equations give relations among generators of the 
path algebra. It is natural to consider the ideal T generated by the F-term equations and define the 
factor algebra A ~ CQ jT. The bifundamental fields obeying the F-term constraints then generate 
a representation of this factor algebra. Namely, classical configurations of the quiver gauge theory 
obeying the F-term constraints are in one-to-one correspondence with A-modules. 

As an example, the algebra A for the conifold geometry contains an idempotent ring C[ei,e2] 
generated by two elements and is given by the following four generators and relationsjf] 

A = C[ei, e2](ai, 02, 61, 62)/ {a\hia2 = a2&iai, 610462 = ^2ado\)i=x.2 ' (3-2) 

Each A-module for this algebra corresponds to a choice of ranks of the gauge groups and a config- 
uration of the bifundamental fields ai^hi satisfying the F-term constraints. 

F-term constraints have a nice geometric interpretation on the quiver diagram, which we will 
find useful in the next section. We observe that each bifundamental field appears exactly twice in 
the superpotential with different signs of coefficients in the superpotential shown in (j2.2p . By taking 
a derivative of the superpotential with respect to a bifundamental field corresponding to a given 
arrow, the resulting F-term constraint states that the product of bifundamental fields around a face 
of the quiver on on one side of the arrow is equal to that around the face on the other side. See 
Figure m for an example. Therefore, when we have a product of bifundamental fields along a path, 
any loop on the path can be moved along the path and the resulting product is F-term equivalent 
to the original one. In [10], it is shown that for any point i, j G Qo, we can find a shortest path Vi^j 
from i to j such that any other path a from i to j is F-term equivalent to Vijuj" with non-negative 

^The center Z{A) of this algebra A is generated by Xij = aihj + bjai(i,j = 1, 2), and is given by 

Z(A) = <C[x-i_i,X-i2,X21,X22\/ {xi-i_X22 - X12X21), (3.1) 
which is the ring of functions of the conifold singularity. 
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integer n, where w is a loop around one face of the quiver diagram. It does not matter where the 
loop Lo is inserted along the path Vij since different insertions are all F-term equivalent. This means 
that any path is characterized by the integer n and the shortest path Vi^j . 



Figure 4: Representation of F-term constraints on the quiver diagram on T^. In this example, if we 
write by Xab the bifundamental corresponding to an arrow starting from vertex A and ending at B 
etc., then the superpotential ET^ contains a term W = —iv{XABXBcXcA)+^T:{XABXBDXDEXEA)^ 
and the F-term condition for Xab (multiplied by Xab) says that the product of bifundamentals 
fields along the triangle ABC and that along the square ABDE is the same. 

In the next subsection, we will impose the D-term constraints on the space of finitely generated 
left A-modules, mod A. Before doing this, however, it is instructive to discuss topological aspects 
of mod A by considering its bounded derived categorj^ D'^ {mod A) . In mathematics, the algebra A 
gives the so-called "non-commutative crepant resolution" [45]. For singular Calabi-Yau manifolds 
such as Xa, the crepant resolution means a resolution that preserves the Calabi-Yau condition|f| 
For a crepant resolution Ya of Xa, we have the following equivalence of categoriej^: 



where D^(coh(YA)) is a bounded derived category of coherent sheaves of crepant resolution Ya, and 
D^{modA) is the bounded derived categories of finitely generated left A-modules. The equation 
()3.3p is also interesting from the physics viewpoint. Since _D''(coh(YA)) gives a topological classifi- 
cation of A branes on the resolved space Ya , the equivalence means that (mod A) also classifies 
D branes, which is consistent with our interpretation above that A-modules are in one-to-one cor- 
respondence with a configuration of bifundamental fields obeying the F-term constraints. 

We should note that the paper [10], which computes the non-commutative Donaldson-Thomas 
invariants for general toric Calabi-Yau manifolds, requires a set of conditions on brane tilings, 
namely on the superpotential. We find that most of their conditions (specified in lemma 3.5 and 
conditions 4.12) are automatically satisfied for any quiver gauge theories for D branes on general 
toric Calabi-Yau manifold. We have not been able to prove that the condition 5.3 also holds in 
general, but it is satisfied in all the examples we know. 

^See [44] for an introductory explanation of derived categories in the context of string theory. 

^Mathematically, we mean a resolution / : X/^ such that ijJy = f*^Xy where ujx and (jJy are canonical 

bundles of X and Y , respectively. For the class of toric Calabi-Yau threefolds, the existence of crepant resolution is 
known and different crepant resolutions related by flops are equivalent in derived categories [46]. 

^This is well-known in the case of the conifold (cf. [47]). For general toric Calabi-Yau threefolds, this is not yet 
mathematically proven as far as the authors are aware of, although there are proofs in several examples [14, 48, 49] . 




(coh(rA)) = £'''(modA), 



(3.3) 
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3.3 D-term Constraints and ^-Stability 

We saw that the derived category D''{modA) of A- modules gives the topological classification of D 
branes in the toric Calabi-Yau manifold Xj\. To understand the moduli space of D branes, however, 
we also need to understand implications of the D-term constraints. This is where the 0-stability 
comes in0 Let 9 £ N^" be a vector whose components are real numbers. Consider an A-module 
M, and recall that this M is decomposed as M — (Bi^Q^Mi with Mi = CiM . The module M is 
called 6'-stable if 

^ 61, (dim CiM') > 0. (3.4) 

iGQo 

for every submodule M' of M ^ When > is replaced by >, the module M is called 6'-semistable. 

In the language of gauge theory, the stability condition (|3.4p is required by the D-term conditions. 
Some readers might wonder why the D-term conditions, which are equality relations, can be replaced 
by an inequality as p.4p . In fact, the similar story goes for the Hermitian Yang-Mills equations. 
There instead of solving the Donaldson- Uhlcnbeck-Yau equations, we can consider holomorphic 
vector bundles with a suitable stability condition, the so-called /i-stability or Mumford-Takemoto 
stability [52,53]. As we mentioned at the end of section 2, it is known that a configuration of 
bifundamental fields is mapped to a solution to the D-term equations by a complexified gauge 
transformation Gc if and only if the configuration is ^-stable. Since each A-module M gives a 
representation of Gc = IlieQo C!L{Ni,C), where GL{Ni,C) is represented by Ali = CiM at each 
node, each A-module specifies a particular Gc orbit. Thus, finding a ^-stable module is the same 
as solving the D-term conditions. 

Up to this point we have not specified the value of 6. Physically, 9's correspond to the FI 
parameters, which are needed to write down D-term equations [50]. Although the Euler number of 
the space of 0-stable A-modules does not change under infinitesimal deformation of 6, it does change 
along the walls of marginal stability [13,14]. The noncommutative Donaldson-Thomas invariant 
defined by [9] is in a particular chamber in the space of 0's. Following [10], we hereafter take 
= (0, 0, . . . , 0). We will comment more about this issue in the final section. 

3.4 D6 Brane and Compactification of Moduli Space 

We have found that solutions to the F-term and D-term conditions in the quiver gauge theory are 
identified with ^-stable A-modules. We want to understand the moduli space of such modules and 
compute its Euler number. 



'^'^The S-stability is a special limit of Il-stability as discussed in [50, 51] . 
'In some literature, an additional condition J^igQo ^ 
satisfied for the choice S = (0, 0, . . . , 0) we choose below 



'^'In some literature, an additional condition J^igQo Si{dim Mi) = is imposed for a choice of 9. This is trivially 
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Since D brane charges correspond to the ranks of the gauge groups, we consider moduh space 
of 6'-stable modules with dimension dim Mi — Ni {i E Qo), which we denote by A4'^{A). To 
compute its Euler number, we need to address the fact that the moduh space of stable ^-modules 
is not always compact. In mathematics literature, the necessary compactification is performed by 
enlarging the quiver diagram by adding one more node in the following way. 

Let us fix an arbitrary vertex i^, and define a new quiver Q — {Qo, Qi) by 

Qo = QoU{*}, Qi = Qi U{a* : * ^ zo}. (3.5) 

Namely, we have added one new vertex * and one arrow * — s- to obtain the extended quiver 
diagram Q. As in the previous case for Q, we can define the path algebra <CQ, the ideal J- generated 
in CQ by and the factor algebra A = CQ/T. Define 9 e by ^ (6*, 1) and define ^-stable 

and semistable A-modules using stability parameter 9. It is shown in lemma 2.3 of [10] that 9- 
semistable A-modules are always 6'-stable, and the moduli space Aif^{A) of 6'-stable modules with 
specified dimension vector N e N'^n+i is compact. 

Adding the extra-node allows us to compactify the moduli space. In the language of D branes, 
this corresponds to adding a single D6 brane filling the entire Calabi-Yau manifold, which is neces- 
sary to interpret the whole system as a six-dimensional C/(l) gauge theory related to the Donaldson- 
Thomas theory. As we mentioned in section 2, the D6 brane serves as a fiavor brane and adds an 
extra node exactly in the way described in the above paragraph. Note that, in the above paragraph, 
the ideal J- is generated by the original ideal J-. In the quiver gauge theory, this corresponds to the 
fact that the fiavor brane does not introduce a new gauge invariant operator to modify the superpo- 
tential. In this way, we arrive at the definition of non-commutative Donaldson-Thomas invariant as 
the Euler number x{-Mf^{A)) of cohomologies. With our identification of A^^(A) with the moduli 
space of solutions to the F-term and D-term conditions, computes the Witten index 

of bound states of DO and D2 branes bound on a single D6 brane (ignoring the trivial degrees of 
freedom corresponding the center of mass of D branes in M^^'^). 

We have chosen a specific vertex io to define the non-commutative Donaldson- Thomas invariant. 
The io dependence drops out in simple cases such as C"^ and conifold, but in general x(-^i^(^)) 
depends on the choice of the io- We note that the quiver gauge theory discussed in section 2 also 
has an apparent dependence on io- There io corresponds to the Qo-type domain where the D6 
brane is located after the T-duality. Since the fundamental group of the toric Calabi-Yau manifold 
Xa is trivial, there is no moduli intrinsic to the D6 brane before taking the T-duality. Thus, we 



expect that the apparent io dependence in the gauge theory side should disa 
treatment of the T-duality. It would be interesting to study this point further 



roear with a proper 
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'^^In the example discussed in [14], the result is dependent on the choice of io and a particular vertex should be 
chosen in order to match with the category of perverse coherent sheaves. We thank Yukinobu Toda for discussions 
on this point. 
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4 Crystal Melting 



In this section, we define a statistical mechanical model of crystal melting and show that the 
model reproduces the counting of BPS bound state of D branes. Using the quiver diagram and the 
superpotential of the gauge theory, we define a natural crystalline structure in three dimensions. 
We specify a rule to remove atoms from the crystal and show that each molten crystal corresponds 
to a particular BPS bound states of D branes. We use the result of [10] to show that all the relevant 
BPS states are counted in this way. 

4.1 Crystalline Structure 

Mathematically, the three-dimensional crystal we define here is equivalent to a set of basis for Aetg , 
where A is the factor algebra A = CQ/J^ of the path algebra CQ divided by the ideal generated 
by the F-term constraints and e,;, is the path of zero length at the reference node io, which is also 
the projection operator to the space of paths starting at zq. Colloquially, the crystal is a set of 
paths starting at io modulo the F-term relations. As we shall see later, it corresponds to a BPS 
state corresponding to a single D6 brane with no DO and D2 charges. We interpret Aeio '^^ terms 
of a three-dimensional crystal as follows. 

The crystal is composed of atoms piled up on nodes in the universal covering Q on M?. By 
using the projection, w : Q ^ Q, each atom is assigned with a color corresponding to the node in 
the original quiver diagram Q. The arrows of the quiver diagram determines the chemical bond 
between atoms. We start by putting one atom on the top of the reference node iq. Next attach an 
atom at an adjacent node j G Qo that is connected io by an arrow going from io to j. The atoms 
at such nodes are placed lower then the atom at zq. In the next step, start with the atoms we just 
placed and follow arrows emanating from them to attach more atoms at the heads of the arrows. 

As we repeat this procedure, we may return back to a node where an atom is already placed. In 
such a case, we use the following rule. As we explained in section 3.2, modulo F-term constraints, 
any oriented path a starting at io and ending at j can be expressed as Vigjoj"', where u) is the loop 
around a face in the quiver diagram and Vi^j is one of the shortest paths from io to j. This defines 
an integer h{a) = n for each path a. The rule of placing atoms is that, if a path a takes io to j 
and if h{a) = n, we place an atom at the n-th place under the first atom on the node j. If there is 
already an atom at the n-th place, we do not place a new atom. 

By repeating this procedure, we continue to attach atoms and construct a pyramid consisting of 
infinitely many atoms. Since atoms are placed following paths from io modulo the F-term relations, 
it is clear that atoms in the crystal are in one-to-one correspondence with basis elements of Aetg . 
Note that by construction the crystal has a single peak at the reference node io- 

This defines a crystalline structure for an arbitrary toric Calabi-Yau manifold. In particular, it 
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reproduces the crystal for C'^ discussed in [2,5], and the one for conifold in [9]. See Figure [5] for the 
crystaUine structure corresponding to the Suspended Pinched Point singularity. In this example, 
the ridge of the crystal (shown as blue lines in Figure [5]) coincides with the {p, q)-weh of the toric 
geometry. As we will discuss later, this is a general property of our crystal. 




Figure 5: Starting from the universal cover Q of quiver Q shown on the left, we can construct a 
crystal on the right. Each atom carries a color corresponding to a node in Q, and they are connected 
by arrows in Qi. The green arrows represent arrows on the surfaces of the crystal, whereas the red 
ones are not. In the case of the Suspended Pinched Point singularity, the atoms come with 3 colors 
(white, black and gray), corresponding to the 3 nodes of the original quiver diagram Q on shown 
in Figure [T] 

4.2 BPS State and Molten Crystal 

In the forthcoming discussions, the crystal defined above will be identified with a single D6 brane 
with no DO and D2 charges. Bound states with non-zero DO and D2 charges are obtained by 
removing atoms following the rule specified below. 

In [9, 10], the Donaldson-Thomas invariants x(7Vl^(^)) are computed by using the [/(l)**^ 
symmetry of the moduli space Mfj^ corresponding to translational invariance of T'^. By the standard 
localization techniques, the Euler number can be evaluated at the fixed point set of the moduli space 
under the symmetry. Correspondingly, in the gauge theory side, BPS states counted by the index 
are those that are invariant under the global [/(l)**^ symmetry acting on bifundamcntal fields 
preserving the F-term constraints since those do not have extra zero modes and do not contribute 
to the index. We are interested in counting such BPS states. 

In order for a molten crystal to correspond to [/(l)®^ invariant ^-stable A- modules, we need to 
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impose the following rule to remove atoms from the crystal. Let f2 be a finite set of atoms to be 
removed from the crystal. 

The Melting Rule: If aa is in for some a € A, then a should also be in fJ. 

Since atoms of the crystal correspond to elements of Aci^ , we used the natural action of A on Aei^ 
to define aa in the above. This means that crystal melting starts at the peak at zq and takes place 
following paths in Aci^ . An example of a molten crystal satisfying this condition is shown in Figure 




Figure 6: Example of a molten crystal and its complement Vt. In this example ^ contains 12 atoms, 
one hidden behind an atom on the reference point represented by a blue point. It is easy to check 
that n satisfies the melting rule mentioned in the text. 

The melting rule means that a complement I of the vector space spanned by in Ae^o gives an 
ideal of A. To see this, we just need to take the contraposition of the melting rule. It states: For 
any /3 e X and for any a G A, a/3 is also in I. 

Generally speaking, an ideal of an algebra defines a module. To see this, consider a vector \I) 
which is annihilated by all elements of the ideal I. From \I) , we can generate a finite dimensional 
representation of the algebra A by acting elements of A on it. However, the converse is not always 
true. Fortunately, when modules are 0-stable and invariant under the U{1)®'^ symmetry, it was 
shown in [9, 10] that there is a one-to-one correspondence between ideals and modules. It follows 
that our molten crystal configurations are also in one-to-one correspondence with ^-modules and 
therefore with relevant BPS bound states of D branes. This proves that the statistical model of 
crystal melting computes the index of D brane bound states. 
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It would be instructive to understand explicitly how each molten crystal configuration corre- 
sponds to a BPS bound state. Starting from a molten crystal specified by f2, Prepare a one- 
dimensional vector space Va with basis vector Bq, for each atom a £ Q. For each arrow a of Q, 
define the action of a on Va by a{ea) = when the arrow a begins from a and ends an another 
atom /3 e J7. Otherwise a{ea) is defined to be zero. Since an arbitrary path is generated by con- 
catenation of arrows, we have defined an action of a £ A on each Va- By linearly extending the 
action of a onto the total space M = (BaenVa, we obtain a A-module M. 

There are several special properties about this module M. First, the F-term relations are 
automatically satisfied. This is because when there exists two different paths a,b £ A starting 
at a and ending at /3, a{ea) and b{ea) are both defined to be 6/3. Second, by construction M is 
generated by action of the algebra A on a single element e^^o £ Vi,o- In such a case AI is called a 
cyclic A-module, and by lemma 2.3 of [10] is also 6'-stable. Third, by the cyclicity of the module it 
follows that M is J7(l)®^ invariant up to gauge transformations. Therefore, M is a J7(l)®^ invariant 
^-stable module. It follows from the result of [10] discussed at the beginning of this section that M 
indeed corresponds a bound state of D branes contributing to the Witten index. 

At the beginning of this subsection, we stated without explanation that the original crystal 
corresponds to a single D6 brane with no DO and D2 charges, and removing atoms correspond to 
adding the D brane charges. To understand this statement, let us recall that, in section 2, we 
started with a configuration of DO and D2 branes on the toric Calabi-Yau manifold and took a 
T-duality along the fiber to arrive at the brane configuration. Thus, the number of D2 branes at 
each node j of the quiver diagram Q is a combination of DO and D2 charges before T-duality. It is 
this number that is equal to the rank of the gauge group at j. 

By using the projection vr : Qo ^ Qo, the A-module M is decomposed as M = ®j^Q^Mj as we 
saw in section 3, where 

M,^ Va. (4.1) 

Q^z^^,7^(Q)— j 

In particular, the formula (|4.ip means that the rank of the gauge group Nj = dimAfj at the node j 
is equal to the number of atoms with the color corresponding to the node j that have been removed 
from the crystal. Thus, removing an atom at the node j is equivalent to adding DO and D2 charges 
carried by the node j. It is interesting to note that each atom in the crystal does not correspond 
to a single DO brane or a single D2 brane, but each of them carries a specific combination of 
DO and D2 charges. In the crystal melting picture, fundamental constituents are not DO and D2 
branes but the atoms. This reminds us of the quark model of Gell-Mann and Zweig, where the 
fundamental constituents carry combinations of quantum numbers of hadrons, as opposed to the 
Sakata model, where existing elementary particles such as the proton, neutron and A particle are 
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chosen as fundamental constituents. 

4.3 Observations on the Crystal Melting Model 

We would like to make a few observations on the statistical model of crystal melting that counts 
the number of BPS bound states of D branes. 

We have studied several examples of toric Calabi-Yau manifolds and found that the crystal 
structure in each case matches with the toric diagram. In particular, the ridges of the crystal, when 
projected onto the plane, line up with the (p, q) web of the maximally degenerate toric diagram. 
This phenomenon is discussed in Appendix. There, we also explain the correspondence between 
molten crystal configurations and perfect matching of the bipartite graph introduced in section 2.2. 

So far, we have considered molten crystals that are obtained by removing a few atoms. We may 
call them low temperature configurations. The high temperature behavior of the model, describing 
bound states with large DO and D2 charges, is also interesting. For C^, it was shown in [2, 5] 
that the high temperature limit of the crystal melting model reproduces the geometric shape of 
the mirror manifold. Since the high temperature limit of a general statistical model of random 
perfecting matchings is known to be described by a certain plane algebraic curve [54], it would be 
interesting to understand its relation to the mirror of a general toric Calabi-Yau manifold [55] . 

In the last subsection, we found it useful to describe BPS bound states using ideals of the algebra 
A. In the case when the toric Calabi-Yau manifold is C^, ideals are closely related to the quantization 
of the toric structure as discussed in [2]. The gauge theory for is the dimensional reduction of 
the TV = 4 supersymmetric Yang-Mills theory in four dimensions down to one dimension, and the 
bifundamental fields are three adjoint fields. The F-term and D-term conditions require that they 
all commute with each other. Thus chiral ring is generated by three elements x, y, z which commute 
with each other without any further relation. In this case, any ideal T^r is characterized by the 
three-dimensional Young diagram tt. Locate each box in the 3d Young diagram n by the Cartesian 
coordinates {i,j,k) {i,j,k — 1, 2, 3, ...) of the corner of the box most distant from the origin, and 
define to be a set of the 3d Cartesian coordinates (i, j, k) for boxes in tt. We can then define 
the ideal u)„ of the chiral ring by, 

In = {x'-'y'-'w>'-'\{i,j,k) i fi,}. (4.2) 

In [2], this description was obtained by quantizing the toric geometry by using its canonical Kahler 
form and by identifying x'^~^y^~^w^~^ as states in the Hilbert space. 

This can be generalized to an arbitrary toric Calabi-Yau manifold as follows. One starts 
with the quiver diagram corresponding to Xa and use the brane tiling to identify the F-term 
equations. This gives the chiral ring generated by bifundamental fields obeying the F-term and 
D-term relations. As we saw in section 4.3, each BPS bound state is related to an ideal of the chiral 
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algebra. Wo expect that such ideals arise from quantization of the toric structure. BPS boimd 
states of D branes emerging from the quantization of background geometry is reminiscent of the 
bubbling AdS space of [56] and Mathur's conjecture on black hole microstates [57]. 

5 Summary and Discussion 

In this paper, we established the connection between the counting of BPS bound states of DO and 
D2 branes on a single D6 brane to the non-commutative Donaldson-Thomas theory. We studied 
the moduli space of solutions to the F-term and D-term constraints of the quiver gauge theory 
which arises as the low energy limit of the brane configuration. We found the direct correspondence 
between the gauge theory moduli space and the space of modules of the factor algebra of the path 
algebra for the quiver diagram quotiented by its ideal related to the F-term constraints, subject 
to a stability condition to enforce the D-term constraints. Using this correspondence, we found a 
new description of BPS bound states of the D branes in terms of the statistical model of crystal 
melting. The crystalline structure is determined by the quiver diagram and the brane tiling which 
characterize the low energy effective theory of D branes. The crystal is composed of atoms of 
different colors, each of which corresponds to a node of the quiver diagram, and the chemical bond 
is dictated by the arrows of the quiver diagram. BPS states are constructed by removing atoms 
from the crystal. 

The relation between the commutative and non-commutative Donaldson- Thomas invariant has 
been extensively discussed in the recent literature. The degeneracy of D brane bound states changes 
when the value of 6, used to define the stability condition, jumps along the codimension one sub- 
space, which is called walls of marginal stability. The jump in the degeneracy can be computed 
by the wall crossing formula [58-60], and if wc start from a particular chamber and applying wall 
crossing formula, we can obtain the value of x{-Mf{A)) in any chamber we want. In the example of 
conifold [13] , wall crossing relates non-commutative Donaldson- Thomas invariants to commutative 
Donaldson- Thomas invariants and to new invariants defined by Pandharipande and Thomas [61]. 
This story is further generalized by [14] when A has no internal lattice point. When the toric 
diagram contains an internal lattice point, non-commutative Donaldson-Thomas invariant includes 
D4 branes, since H4{Ya) 7^ 0. Since (commutative) Donaldson-Thomas invariants does not include 
D4 brane charges, the above discussion of wall crossing should be modified. 

It has been proven recently that the topological string theory is equivalent to the commutative 
Donaldson- Thomas theory for a general toric Calabi-Yau manifold [2,8]. Since the commutative 
Donaldson- Thomas theory count BPS states for some choice of stability condition, 

Zbh = Ztop, (5.1) 
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is indeed true in some chamber. On the other hand, our result shows that the relation, 




^crystal melting ? 



(5.2) 



holds in another chamber, where Zgjj is the BPS state counting for another choice of the stability 
condition. Combining these two results, we find that the topological string theory and the statistical 
model of crystal melting are related by the wall crossing, and we have 



Since there is no wall crossing phenomenon for the Donaldson-Thomas theory on C'^, this result 
does not contradict with [2], where a direct identification of the topological string theory and the 
crystal melting is made for C"^. In general, we expect that a proper understanding of the relation 
between the topological string theory and the crystal melting requires that we take the wall crossing 
phenomena into account. 

The OSV formula (1.1) suggests yet another relation between the black hole microstate counting 
and the topological string theory. According to [62], for a compact Calabi-Yau manifold, the 
D6/D2/D0 brane system gives rise to a large black hole in four dimensions since it is related to 
a spinning M theory black hole by the Kaluza-Klein reduction. In fact, one can compute the 
semi-classical Bekcnstcin-Hawking entropy for such a 4-dimensional black hole and find that it can 
be made arbitrarily large provided the D2 charges are sufficiently larger than the DO charge. In 
this paper, we discussed the D6/D2/D0 system on a non-compact toric Calabi-Yau manifold with 
an infinite volume. Though it is not obvious that the gravity description in four dimensions is 
applicable in this case, the OSV formula has been successfully tested for a similar class of non- 
compact Calabi-Yau manifolds [63 66]. If it is applicable in our case, it would imply the relation 
between .^crystal melting and l^topp, modulo wall crossings. It would be interesting to find out if 
such a relation holds. 

It appears that the crystal melting picture is closely related to the quantization of the toric 
structure of the Calabi-Yau manifold. It would be interesting to understand the relation better. 
This could lead to a new insight into quantum geometry, along with the observations in [56] for the 
bubbling AdS geometry and [57] for black hole microstates. 
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A Perfect Matchings 

In this Appendix we are going to explain the one-to-one correspondence between a molten crystal 
discussed in the main text and a perfect matching of the bipartite graph. This means that the 
problem of counting BPS states can also be reformulated as a problem of counting perfect matchings 
of the bipartite graph, where a perfect matching is a subset of edges of the bipartite graph such that 
each vertex is contained exactly once. The contents of this appendix is basically a recapitulation 
of [10]. 

In section [4.11 we considered a quiver Q = {Qo, Qi), which is a universal cover of the quiver Q 
on T^. The dual graph of Q, which we denote by P, can be made bipartite using orientation of 
arrows of Q, and is a universal cover of the bipartite graph P on described in section 2.2. What 
we are going to do is to give an explicit correspondence between a perfect matching of the bipartite 
graph P and a configuration of molten crystal. 

We first construct a perfect matching from a molten crystal. Given a molten crystal as shown 
in Figure O choose all the arrows of Q which are along the surface of the crystal. In the example 
of Figure [5l such arrows are colored green in Figure [5l while the remaining arrows are colored red. 
Take the set of the dual of edges colored red. It is proven by [10] that such a subset of edges of P 
is a perfect matching. This is the perfect matching we wanted to construct. 

In the case when no atoms are removed from the crystal, the perfect matching obtained by this 
method is called the canonical perfect matching, which we denote by Dq. Since only a finite number 
of atoms are removed from the crystal, the perfect matching obtained from a molten crystal by the 
above method coincides with D when sufficiently away from the reference point i. 

Conversely, given a perfect matching D which coincides with Dq when sufficiently away from 
the reference point i, we can reproduce a molten crystal. Let us superimpose D with Dq, and we 
have a finite number of loops, as shown in Figure [5] in the case of Suspended Pinched Point. Define 
a height function hi) such that 

(1) h]j{j) ~ when sufficiently away from i. 

(2) h]j increases by one whenever we cross the loop and go inside it. 

The example of ho for the case of Suspended Pinched Point is shown in Figure [51 

By removing h]j{j) atoms from each j e Qo, we can construct a molten crystal. It was proven 
in [10] that the set of atoms removed from the crystal so defined satisfies the melting rule of section 



19 



Figure 7: Given a configuration of a molten crystal, we can construct a perfect matching of the 
bipartite graph. Each arrow is colored green if it is along the surface of the crystal, and red 
otherwise. The set of dual of arrows colored red gives a perfect matching of the bipartite graph. 




Figure 8: By superimposing a perfect matching of Figure [7] with the canonical perfect matching 
shown later in Figure [H we have a set of loops, which defines a height fimction Hd. From this 
function we can recover a molten crystal. 
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14.21 This establishes the one-to-one correspondence between a molten crystal and a perfect matching 
of the bipartite graph, meaning that BPS states can also be counted by perfect matchings of the 
bipartite graph T. 

Finally, let us finish this Appendix by pointing out an interesting connection of the canonical 
perfect matching Dq with toric geometry. The example of canonical perfect matching Dq for the 
Suspended Pinched Point is shown in Figure[51 In this example, the asymptotic form of the bipartite 
graph has four different patterns. Each of four patterns is periodic and therefore an be thought 
of as a perfect matching of the bipartite graph on T^. In the brane tiling literature, a perfect 
matching on the bipartite graph on is known to correspond to one of the lattice points of the 
toric diagram [25,4O]0 We recognize that the four perfect matchings are identified with the four 
corners of the toric diagram in Figure l-(a) and that the borders between different patters are 
identified with the blue lines in Figure l-(b), which makes the (p, q)-'weh of the diagram. 

In general, for an arbitrary toric Calabi-Yau manifold, we can use the same pattern to construct 
a perfect matching. Divide the universal covering of the bipartite graph into segments separated 
by the (p, q)-'weh of the toric diagram^ The perfect matching in each segment is periodic and is 
identified with one of the perfect matchings of bipartite graphs on T^, which corresponds to one 
of the lattice points of the toric diagram and the lattice point in question is precisely the vertex 
surrounded by the two {p, q)-weh8 on ^ This determines a perfect matching. In particular, this 
means that the ridges of the crystal line up with the (p, g)-web of the toric diagram. 

We have examined several other examples as well, and this pattern holds in all cases. Thus, we 
conjecture that the perfect matching constructed in this way is canonical. We would like to stress 
again that this conjecture is not needed to construct the crystal melting model. Here we are simply 
pointing out that, in the examples we have studied, the crystalline structures fit beautifully with 
the corresponding toric geometries. 
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